String effects in Polyakov loop correlators 
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Abstract 

We compare the predictions of the effective string description of confine- 
ment in finite temperature gauge theories to high precision Monte Carlo 
data for the three-dimensional Z2 gauge theory. First we review the pre- 
dictions of the free bosonic string model and their asymptotic behavior 
in the various regimes of physical interest. Then we show that very good 
agreement with the Monte Carlo data is obtained, for temperatures not 
too close to the deconfinement one (typically T < T c /3). For higher tem- 
peratures, higher order effects are not negligible: we show that they are 
accurately modeled by assuming a Nambu-Goto string action and com- 
puting its partition function at next-to-leading order. 

1 Introduction 

Two color sources in a confining gauge theory are bound together by a thin flux 
tube, which can fluctuate like a massless string: this hypothesis constitutes the 
core of the effective string description of confinement. Besides providing an ap- 
pealing, qualitative explanation for the linear behavior of the confining potential 
at large distances, this picture provides us with quantitative, testable predic- 
tions about this highly non perturbative regime of gauge theories, as pointed 
out in the seminal papers by Luscher, Symanzik and Weisz jl], |2j: the quantum 
fluctuations of the flux tube produce measurable effects on the gauge invariant 
correlation functions of the gauge theory, namely Wilson loops and Polyakov 
loop correlation functions. 



Among these effects, the most widely known is the Liischer term correction 
to the zero-temperature linear confining potential: 

V(R) ~ aR - (1) 

where c is the conformal anomaly of the two-dimensional field theory describing 
the flux tube fluctuations: c = d — 2 for a free bosonic string in a d-dimensional 
spacetime. 

The predictive power of the effective string picture goes well beyond Eq. (|l|) , 
as it extends to the full functional form of the Wilson loop expectation values 
at large distances: 

< W(R, T) >= e- aRT+ ^ R+T ^ +k Z q (R, T) (2) 

where Z q (R,T) is the partition function of the two-dimensional quantum field 
theory describing the quantum fluctuations of the flux tube. For the free bosonic 
string, this is simply the theory of d — 2 free massless scalar fields living on the 
rectangle defined by the Wilson loop: 
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where T)(t) is the Dedekind r\ function and r = iT/R. The prediction Eq.(^|) 
was successfully compared to Monte Carlo data for three-dimensional Z2 gauge 
theory in Ref.||, thus confirming the validity of the string description, and 
showing its full predictive power. 

It must be kept in mind that the free string description that gives Eq.(^|) 
is an effective description that is expected to hold in the long distance regime. 
The "true" string theory describing gauge theories at all length scales, if it 
exists at all, certainly contains complicated string self-interactions, probably 
described by non-renormalizable terms in the string action. Precisely because 
they are non-renormalizable, however, these terms should become negligible in 
the infrared limit, so that the free string model is expected to describe the 
physics of confinement for large distances independently of the details of the 
full, interacting theory. Among the results of Ref . fl is a quantitative estimate 
of the typical distances, in physical units, where the free string picture becomes 
numerically accurate: if the physical distances studied are too small, one cannot 
unambiguously determine which string model actually describes the flux tube 
fluctuations^]. 

In this paper we consider the effective string description for a confining gauge 
theory at finite temperature. The gauge invariant quantity of interest is now 
the correlation function of two Polyakov loops. The effective string description 
predicts its dependence on the temperature and on the distance between the 

1 This is the reason why earlier studies, probing shorter physical distances due to the smaller 
computational power available, could not identify the free bosonic string as the correct model, 
and actually suggested a fermionic stringmodel as the best description of the confining string 
both at zero Q and finite temperature jij. 
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two loops. It is this prediction that we want to compare to Monte Carlo data. 
We chose to work again with the three-dimensional Z2 gauge theory, where very 
high precision can be achieved on large lattices, to make our test as stringent 
as possible. 

Also in the finite-temperature case it is important to keep in mind the fact 
that the free string picture is an effective description valid in the infrared limit, 
and to determine quantitatively the distance and temperature scales where it 
becomes numerically accurate. First of all, one has to take into account the 
known fact that the actual color flux tube has a finite thickness of the order 
~ 1/T C , where T c is the deconfinement critical temperature therefore the 
free string description, which is based on an idealized, one-dimensional flux tube, 
certainly must break down at distances lower that 1/T C . Moreover, we will show 
that the string self-interaction effects mentioned above cannot be neglected for 
temperatures close to T c : within the accuracy of our data, one has to take 
T < T c /3 to find good agreement between Monte Carlo data and free string 
predictions. 

We will also show that by modeling the self-interaction of the confining 
string by a Nambu-Goto action and computing the string partition function at 
next-to-leading order one can successfully take into account the most important 
corrections to the free string picture: one then obtains good agreement with the 
Monte Carlo data for temperatures up to ~ T c /2 (without introducing any new 
adjustable parameters). 

Since the effective string description is believed to be universal, that is to 
hold for all confining gauge theories irrespective of the gauge group or the space- 
time dimensionality, our conclusions are likely to apply also to SU(N) gauge 
theories in 2 + 1 and 3 + 1 dimensions. Agreement between free string theory 
predictions and lattice determinations of the interquark potential in SU (N) 
gauge theories was recently reported by various groups R |8|, ^] . 

2 Effective string predictions 

The quantum fluctuations of the flux tube at finite temperature are described, 
in the free string model, by d — 2 free massless scalar fields obeying Dirichlet 
boundary conditions on the two Polyakov loops and periodic boundary con- 
ditions in the compactified direction (euclidean time). The correlation of two 
Polyakov loops at distance R is then predicted to be 



where the free energy F(R,L) depends on the inverse temperature L = 1/T 
(i.e. the lattice size in the time direction) and the distance R, and is given by 
a classical and a quantum contribution: 



(P(O)pt(fl)) = exp l-F(R,L)} 



(4) 



F(R,L) = F cl (R,L) + F q (R,L) 
The classical term corresponds to the area law: 
F c i(R,L)=a LR + k(L) 



(5) 



(6) 
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Since the difference between Wilson loops and Polyakov loop correlation func- 
tions is, in this picture, completely taken into account by the different choices 
of boundary conditions for the string fluctuations, <7q is the zero-temperature 
string tension, extracted, say, from the Wilson loop at the same value of the 
gauge coupling. 

The quantum term encodes the flux tube fluctuations, and is equal to the 
free energy of the d — 2 massless scalar fields describing them|T(| [II], |l2| : 



F q (R,L) = {d -2) log r)(r) t 

where 77 is again the Dedekind eta function: 

00 

V{r) = Y[(l - q n ) ] q = e 
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According to the value of the ratio — it = L/2R, one can use the two expansions: 
2R < L 



F q (R,L) = (d-2) 
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The latter expression shows that the Polyakov loop correlation function should 
decay at large R as 

a(T) 



(P(0)P'(R)) ~ exp ; —R 

with a temperature-dependent string tension a(T) given by 

ttT 2 



a(T) = a -(d-2)- 



G 



(11) 
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3 Comparison with Monte Carlo data 

3.1 At what scales is the free string picture expected to 
hold? 

To compare the free string predictions described in the previous section to Monte 
Carlo data, we need to establish in which regime of distance R and temperature 
1/L we expect these predictions to be fulfilled. 

It is known that the picture of a one-dimensional confining string is an 
idealization valid for long distances. The flux tube has actually a finite thickness 
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of order 1/T C , where T c is the deconfinement critical temperature Q. Therefore 
we expect the free string picture to be accurate for R > 1/T C . 

We also know that the free string picture must break down for temperatures 
close to T c : in fact if we were to believe the free string picture for all tempera- 



tures up to T c , Eq.([12[) would predict the value of the latter to be T c = y/6<7o/n, 
a prediction that turns out to be very far from the true value. 

A much better prediction for the dimensionless ratio T^/ao is obtained if 
one assumes for the confining string a Nambu-Goto action, proportional to the 
area of the surface spanned by the string. One can then predict || 



a(T) = a o yi-^-j (13) 
and 

which reduces to Eq.([l2]) for T « T c . The impressive agreement one finds 
with Monte Carlo data for various gauge groups suggests that these equations 
give at least a good approximation to the true behavior. 

Therefore we can use Eq. (|l3| ) to estimate the range of temperatures in which 
string interaction effects can be safely neglected and the free string picture is 
expected to be accurate: expanding the square root in Eq.([l3|) to next-to- leading 
order we find 



a(T) ~ a 



1 (T\ 1 (T 



2\T C J 8 \T. 



(15) 



We expect the free string picture to give an accurate description of the data when 
the last term in parentheses is comparable to the accuracy in the determination 
of (To. In our case such accuracy is of order 1%, so that we do not expect to 
be able to use the free string prediction Eq.(|l2|) for temperatures higher than 
~ T c /2. We will see below that this estimate turns out to be too optimistic, since 
good agreement between free string and Monte Carlo data is obtained only up to 
~ T c /3. Note however that this determination of the limiting temperature ratio 
T/T c for the free string picture depends in an essential way from the accuracy 
of the Monte Carlo data, and hence has no intrinsic meaning. 

3.2 Monte Carlo simulation 

We simulated the three-dimensional Z2 gauge model with Wilson action: 

S = - ^ <7 l f7 20'3Cr4 (16) 

□ 

where the sum is extended to all plaquettes of a cubic lattice, and a%, . . . ,(74 
are Z2 variables defined on the four links around the plaquette. The partition 
function is 

Z(f3) =^exp-/3S (17) 
M 
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where the coupling (3 is the same for all directions. 

With these conventions the model, at zero temperature, is known to undergo 
a roughening transition at (3 r = 0.4745(2), where the strong coupling expansion 
of the Wilson loop ceases to converge, signaling the fact that the flux tube 
fluctuations become massless: for (3 > (3 r we expect the free bosonic string 
description to hold. Increasing (3 we find a (bulk) deconfinement transition at 
f3 c = 0.76141337(22) where the string tension extracted from Wilson loops 
vanishes. 

At finite temperature, the model undergoes a deconfinement transition, at 
a coupling (3d{L) that depend on the inverse temperature L = l/T. For several 
values of L the critical coupling is known to high precision. In particular we 
have 

(3 d {A) = 0.73107(2) (18) 
/3 d (6) = 0.746035(8) (19) 

In our simulations we used the gauge version of the microcanonical demon 
algorithm Jl1| [hJ combined with the canonical update of the demons ||(| , de- 
scribed in detail in Ref. pifl , and implemented, as in Ref. §3], in multi-spin 
coding technique. We fixed the coupling (3 to coincide with one of the values 
reported in Eqs. ( |l8||l9| ), so as to be able to control the value of T/T c by vary- 
ing the lattice size L. The lattice size in the two space-like directions was fixed 
at jV s = 48, corresponding to at least ten times the correlation length in the 
space-like directions. 

For these two values of (3 we simulated the model at several values of the 
inverse temperature L, all corresponding to temperatures up to T c /2. In Tab. 2 
we show for each of the two (3 values the values of L we used in the simulations, 
and the value of the zero-temperature string tension <7o(/3). The latter were 
evaluated by interpolating the high-precision estimates published in Refs.p2|, 
p3| p4~f for several reference values of (3 with the scaling formula 

\Pref ~ (3c J 

where f3 re f is the closest value for which a direct Monte Carlo estimate is avail- 
able, and v is the correlation length critical index. Neglecting corrections to 
scaling introduces a systematic error that can be estimated by evaluating o~o {J3) 
using two nearby reference values. The errors quoted in Tab. 2 are the sum of 
the statistical and systematic errors. 

For each (3 and L in Tab. 2 we evaluated the ratio 

where G(R) is the Polyakov loop correlation function at distance R: 

G{R) = {P{x)P^{x + R)) (22) 
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1/T C 


L 


c a 2 


0.73107 


4 


8-14 


0.0440(3) 


0.74603 


6 


12,14,15,16,17,18,20 


0.018943(32) 



Table 1: Temporal lattice sizes and zero-temperature string tension for our 
values of (3. 



for all R between and N s /2 = 24. Errors on Q(R) were evaluated using a 
standard jackknife procedure. 

From Eq. (^3j) we see that if the quantum fluctuations of the confining string 
are neglected, Q{R) is predicted to be a constant: 

Q(R) = a o (0)L (23) 

so that the difference 

Q q {R) = Q(R) - a ([3)L (24) 

can be interpreted as the quantum contribution: the free string prediction for 
this quantity is 

Q q (R) = F q (R+l,L)-F q (R,L) (25) 

with F q given by Eq. (0) . 

The complete Monte Caro results will be presented in a forthcoming publi- 
cation, together with details about the data analysis procedure. Here we limit 
ourselves to a smaller sample of data, sufficient to evidentiate the most impor- 
tant results of the analysis, that are likely to be relevant also for SU(N) gauge 
theories. 



3.3 Results at T < T c /3: the free string. 

All the data at T < T c /3 are in good agreement with the free string predictions: 
in Figs. 1-2 we have plotted the MC values of Q q (R) defined by Eq.(p4|) and 
the free string prediction Eq. ( p5| ) for our two (3 values at T — T c /3. The data 
are plotted as a function of 

i 2R 

, = - = -, (26) 

which is the "natural" variable for the effective string prediction. In the three 
figures we have shown only the data with R > 1/T C) since for lower R we do 
not expect the free string to describe the data. Moreover we have not shown 
the data at very large R, where the statistical error on Q q becomes larger than 
the effect we are studying. 

A few interesting things can be learned by looking at these figures: 

• The data we can use for the comparison are all in the region z ~ 1 . This 
is unavoidable since small values of z correspond to distances smaller than 
the width of the flux tube, as explained above, while data at large z have 
large statistical uncertainties. 
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Figure 1: Comparison between free string prediction and Monte Carlo data for 
Q q {R) as defined in Eq.(||). The Monte Carlo data are taken at /3 = 0.73107 
and T — T c /3 = (12a) , and are plotted as a function of z = 2R/L. The solid 
line is the prediction Eq. ( p^ j for this quantity. The dashed line is the asymptotic 
value — 7r/6L of the free string prediction for large z. 
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Figure 2: Same as Fig. 1 for (3 = 0.74603 and T = T c /3 = (18a)" 1 
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• In this region the free string prediction is very far from its asymptotic 
behavior for large R, where it approaches the constant 

F q (R + 1, L) - F q (R, L)^- {d ~ £ 2)7r (27) 

shown for comparison in the figures: the string correction for z ~ 1 has 
actually the opposite sign with respect to its asymptotic behavior. 

• The function F q (z = 2R/L) has a zero in z ~ 1.8. This is due to the 
fact that the asymptotic term ir/(6L) and the subdominant logarithmic 
term that appears in Eq.(|i"o|) have opposite signs. This in turns makes the 
string corrections much harder to detect in the Polyakov loop correlation 
functions than in the Wilson loop case, where no such cancellation occurs. 



3.4 Results at T ~ T c /2: the interacting string. 

If the same comparison between Monte Carlo data and free string predictions is 
performed at temperatures closer to the deconfinement one, significant discrep- 
ancies begin to appear, showing that the free string contribution is not sufficient 
to account for the quantum fluctuations of the flux tube at such temperatures: 
string interaction effects, or more precisely the self-interaction of the world-sheet 
fields describing the string configuration, become non negligible. 

The simplest way to estimate the effect of these self-interactions is to assume 
a Nambu-Goto action for the effective string. On one hand, this is exactly the 
assumption that leads one to the remarkably successful predictions Eqs. ([T^ , [l^ ). 
On the other hand this same choice of string action accurately describes the 
string interaction effects for fluctuating interfaces in three-dimensional statisti- 
cal models]^, [2E|| . 

From this assumption, one can compute the string partition function to next- 
to-leading order in the dimensionless expansion parameter (aRL)^ 1 , with the 
choice of boundary conditions relevant to the Polyakov loop correlation function, 
namely Dirichlet on the two loops and periodic in the Euclidean time direction. 



This calculation was performed in Ref. 26 using £- function regularization; it 
can be shown |^7j that the same result is obtained with several other choices of 
regularization. 

Including these contribution the string free energy Eq.(Q) becomes 

1 \ 



F^ NLO \R,L) = (d-2) 



^SV(r) - -^-^ [2E 4 (r) - EUr)] 



-O 



(aLR) 2 



1152 aR 3 

(28) 

where Ei and E± are the Eisenstein functions. The latter can be expressed in 
power series: 



E 2 (t) = 1-24^^(71)?" (29) 
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Figure 3: Data for Q q at (3 = 0.73107 and T = T c /2 = {8a)- 1 . The solid line 
is the full NLO string prediction, while the dash-dotted line is the free string 
prediction. 



E 4 {t) = l + 240^a 3 (n)g" 



(30) 



q = e— (31) 

where <r(n) and (73(71) are, respectively, the sum of all divisors of n (including 1 
and n), and the sum of their cubes. The modular transformation properties 

/ 1 \ it / 

(32) 



(33) 



are also useful. Note that the inclusion of next-to-leading terms does not require 
the introduction of any new free parameter, so that the predictive power is the 
same as for the free string case. 

In Figs. 3-4 we show the Monte Carlo data for the string fluctuation con- 
tribution defined in Eq.(24) for T = T c /2 compared to both the free string 
prediction and the Nambu-Goto string at next-to-leading order Eq.(p8|). We 
can see that Eq.(p8|) accurately describes the deviations of the data from the 
free string picture. 

A word of caution is however in order: while the free string prediction Eq.(^) 
is the universal infrared limit of a large class of interacting string models (see 
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Figure 4: Same as Fig. 4 for (3 = 0.74603 and T = T c /2 = (12a)- 1 

e.g. the various models considered in p(|), the interaction contributions depend 
on the choice of the model. Therefore while the validity of the free string model 
in the infrared limit R, L — > oo is basically a consequence of the masslessness 
and bosonic character of the fields describing the string fluctuations, the choice 
of the Nambu-Goto action is a much stronger statement, that cannot as far as 
we know be based on simple physical grounds. Nevertheless, it is certainly note- 
worthy that the same choice of action gives accurate predictions for fluctuating 
interfaces [^3[ and Polyakov loop correlation functions. 

4 Conclusions 

The study of Polyakov loop correlation functions in the three-dimensional Z2 
model allows us to gain some insight on the effect of the fluctuations of the 
confining flux tube, which is likely to be relevant also in more complicated and 
realistic gauge theories. The main points we have elucidated are: 

1. The quantum fluctuations of the confining string give a computable and 
numerically relevant contribution to the Polyakov loop correlation func- 
tions in the confined phase. 

2. The free bosonic string model describes well such contributions provided 
the temperature T is not too close to the critical one: in our case T < T c /3 
is required to find good agreement. 
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3. When the temperature gets closer to T c , next-to-leading effects in the 
string partition function become important. Assuming a Nambu-Goto 
action for the confining string one can compute these effects and find 
excellent agreement with the Monte Carlo data also at T ~ T c /2. 
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